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The crossover between two customary limits of phonon-assisted tunneling, the adiabatic and 
antiadiabatic regimes, is studied systematically in the framework of a minimal model for molecular 
devices: a resonant level coupled by displacement to a localized vibrational mode. Conventionally 
associated with the limits where the phonon frequency is either sufficiently small or sufficiently 
large as compared to the bare electronic hopping rate, we show that the crossover between the two 
regimes is governed for strong electron-phonon interactions primarily by the polaronic shift rather 
than the phonon frequency. In particular, the perturbative adiabatic limit is approached only as 
the bare hopping rate F exceeds the polaronic shift, leaving an extended window of couplings where 
r well exceeds the phonon frequency and yet the physics is basically that of the antiadiabatic 
regime. We term this intermediate regime the extended antiadiabatic regime. The effective low- 
energy Hamiltonian in the traditional and the extended antiadiabatic regime is shown to be the 
(purely fermionic) interacting resonant-level model, with parameters that we extract from numerical 
rcnormalization-group calculations. The extended antiadiabatic regime is followed in turn by a true 
crossover region where the polaron gets progressively undressed. In this latter region, the phonon 
configuration strongly deviates from a simple superposition of just one or two coherent states. The 
renormalized tunneling rate, which serves as the low-energy scale in the problem and thus sets the 
width of the tunneling resonance, is found to follow an approximate scaling form on going from the 
adiabatic to the antiadiabatic regime. Charging properties are governed by two distinct mechanisms 
at the extended antiadiabatic and into the crossover region, giving rise to characteristic shoulders in 
the low-temperature conductance as a function of gate voltage. These shoulders serve as a distinct 
experimental fingerprint of phonon-assisted tunneling when the electron-phonon coupling is strong. 

PACS numbers: 71.38.-k,85.65.+h,72.10.Di 
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I. INTRODUCTION 

The promise of molecular electronics has focused enor- 
mous interest on molecular devices X Typically, such de- 
vices consist of an individual molecule trapped between 
two leads, in-between which a voltage bias is applied. By 
measuring the current flowing across the molecular bridge 
one can investigate the molecule's internal degrees of free- 
dom and their coupling to the leads, which can lead in 
turn to complex many-body physics. While lacking the 
exquisite design and control capabilities of semiconduc- 
tor quantum dots, molecular devices can be produced in 
large quantities, thus allowing for many samples to be 
scanned in a relatively short period of time. At the same 
time, it is not always clear if a molecule has been success- 
fully trapped between the leads, or whether it might have 
been damaged or distorted in the course of preparation. 

From a basic-science perspective, single-molecule tran- 
sistors (SMT) offer two major advantages over their semi- 
conductor counterparts. First, the relevant energy scales 
are notably larger in SMTs, rendering these scales more 
accessible to experiments. Second, the electronic degrees 
of freedom are generally coupled to nuclear vibrational 
modes, providing an extraordinary opportunity to study 
the electron-phonon coupling at the nano-scale. The 
same picture applies to suspended carbon nanotubes, 
where the motion of electrons is coupled to vibrations 
of the tube (vibrons). Indeed, phonon-assisted tunneling 
can lead to a plethora of interesting phenomena, includ- 



ing the appearance of inelastic steps and peaks in the dif- 
ferential conductance^^ the Frank-Condon blockade^ 
and the interplay with the Kondo effect £r— While most 
of the experiments cited above are in general accord with 
theoretical expectations, some issues, such as the sign 
of the inelastic steps at integer multiples of the phonon 
frequency, remain under debatei 11 ' 12 Other theoretical 
predictions awaiting experimental verification include un- 
orthodox variants of the Kondo effect ) 13 ' 14 signatures of 
pair tunneling, 15 and interesting nonequilibrium effects 
on the phonon distribution function, 16 to name but a 
few. 

While much of the theoretical activity on molecular 
devices is presently centered on finite-bias transport, in 
this paper we focus on thermal equilibrium and address a 
specific question pertaining to the nature of the crossover 
between two customary limits of phonon-assisted tunnel- 
ing, the so-called adiabatic and antiadiabatic regimes. 
These two terms are broadly used in the context of 
electron-phonon coupling to indicate the limits where the 
bare electronic motion is either sufficiently fast (adia- 
batic limit) or sufficiently slow (antiadiabatic limit) as 
compared to the phonon vibrations. In the more spe- 
cific context of resonant phonon-assisted tunneling, the 
relevant measure of electronic motion is given by the tun- 
neling rate. Hence, the two limits correspond to whether 
the bare electronic tunneling rate T is either sufficiently 
small or sufficiently large as compared to the phonon fre- 
quency U)q. 
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When r <C wo (we work with units in which h = 1), 
the phonon can efficiently respond to hopping events by 
forming a polaron, suppressing thereby the electronic 
tunneling rate. This suppression, which can be quite dra- 
matic, is manifest, e.g., in a narrowing of the tunneling 
resonance. In the opposite limit, loq "C T, the phonon 
is too slow to respond to the frequent tunneling events, 
having little effect on their rate. Each of these extreme 
limits is rather well controlled theoretically, either in the 
framework of the Lang-Firsov transformation^ 7 - or using 
ordinary perturbation theory in the electron-phonon cou- 
pling. Far less understood is the crossover region between 
the two limits, which lacks a small parameter. 

This general picture neglects, however, another impor- 
tant energy scale: the harmonic potential energy asso- 
ciated with the relative displacement of the phonon be- 
tween different molecular electronic configurations. In 
the case of a single spinless level with the dimension- 
less displacement coupling A [see Eq. ([TJ for an explicit 
definition of A], the so-called polaronic shift is given by 
E p = A 2 wo- Thus, depending on the magnitude of A, the 
polaronic shift E p may exceed T and/or luq. Quantum 
mechanically it is natural to associate E p with a new time 
scale r p = 1/Ep, which may either be shorter or longer 
than the electronic dwell time T,i W eii = 1/T and the pe- 
riod of oscillations r osc = 2tt/ujq. Whether r p has the 
status of a true physical time scale is not immediately 
clear. It can not have real significance for A <C 1, when 
the phonon displacement is small as compared to its zero- 
point motion. Neither does r p play any role for a classical 
oscillator, whose period is independent of the amplitude 
of oscillations. At the same time, E p does show up as an 
additional energy scale for phonon- assisted tunnelingjiS 
although its significance, let alone its role in defining the 
physical boundaries of the adiabatic and antiadiabatic 
regimes, have never been quite resolved. 

From this brief discussion it is clear that the most inter- 
esting and yet most challenging regime is that of strong 
electron-phonon coupling, 1 <C A, where r p , whatever its 
physical interpretation might be, can potentially assume 
a dominant role. The interest in strong electron-phonon 
coupling is further amplified by recent reports of large 
values of A (including some in excess of 5) in suspended 
carbon nanotubes. 7 In the antiadiabatic limit, the tun- 
neling electrons experience strong polaronic dressing for 
1 -C A, reflected in an exponential suppression of the 
renormalized tunneling rate from T to r e g = re <T. 
This dramatic effect raises several basic questions: 

1. When tuning the bare tunneling rate from weak 
(r <C ujq) to strong (luq -C T) coupling, which phys- 
ical parameters set the scale for first leaving the 
polaronic physics of the antiadiabatic regime and 
then entering the perturbative physics of the adi- 
abatic regime? Are these two transitions governed 
by a single scale or are there perhaps two distinct 
scales? 

2. How does the ratio r c g/r evolve from its exponen- 



tially small value e in the antiadiabatic regime 
back to approximately one in the adiabatic limit? 
In other terms, how does the polaron get un- 
dressed? 

3. Are there any distinct experimental signatures of 
the crossover regime that can be detected? 



A. Preliminaries 

In this paper, we answer these questions in detail in the 
framework of the resonant-level model with an additional 
displacement coupling to a localized vibrational mode 
[see Eq. (JTJ)]. Besides being one of the most widely used 
models for single-molecule devices, our motivation for 
adopting this specific Hamiltonian is two-fold. The first is 
physical clarity, as this Hamiltonian constitutes the min- 
imal model where the crossover from the antiadiabatic to 
the adiabatic regime of phonon- assisted tunneling can be 
studied without being masked by competing many-body 
effects (e.g., the Kondo effect in case of a spinful level). 
The second point is technical in nature, and pertains to 
our method of choice for accurate nonperturbative cal- 
culations. In this work we employ Wilson's numerical 
renormalization-group (NRG) approach ; 19 i 20 which is a 
highly precise tool for calculating equilibrium properties 
of quantum impurity systems. In the NRG, the com- 
putational effort grows exponentially with the number 
of conduction-electron species, hence the restriction to a 
single spinless band allows us to accurately address large 
values of A that otherwise would be inaccessible using 
more elaborate models. 

We note in passing that the NRG has been applied 
to moderately large values of A in the framework of the 
spinful Anderson-Holstein model^ however, these stud- 
ies used rather large values of the phonon frequency. The 
combination of large values of A and small phonon fre- 
quencies (as is appropriate for molecular and nanotubc 
devices) is presently a challenge to treat accurately using 
the NRG if spin is to be included. 

Focusing on 1 <C A, Fig. Q] displays our resulting sce- 
nario for the crossover from the antiadiabatic to the 
adiabatic regime. In contrast to the naive picture, the 
crossover is governed primarily by the polaronic shift E p 
rather than the phonon frequency ujq. In particular, the 
perturbative adiabatic limit is approached only as T ex- 
ceeds -Ep, leaving a rather broad window of couplings 
where T well exceeds wo and yet the physics is basically 
that of the antiadiabatic regime. We term this interme- 
diate regime the extended antiadiabatic regime. 

The physical origin of the extended antiadiabatic 
regime is in the enormous separation of scales between 
r and T eS when r <C E p . Although smaller than the 
bare tunneling rate, the phonon frequency ujq remains no- 
tably larger than the renormalized one, hence the phonon 
can still efficiently respond to the individual tunneling 
events. The distinction between the extended and the 
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FIG. 1: (Color online) Schematic description of the crossover 
from the antiadiabatic to the adiabatic regime of resonant 
phonon-assisted tunneling, for strong electron-phonon inter- 
actions 1<A. Upon increasing F, an extended antiadiabatic 
regime persists until the renormalized tunneling rate r c fl ap- 
proaches the phonon frequency ujo- The condition r e fi = cjo 
is typically met for T ~ 0.55-E p with E p — \ 2 u>o, extend- 
ing the physics of the antiadiabatic regime far beyond the 
traditional condition V < cjq. In terms of the bare model 
parameters, the extended antiadiabatic regime persists up to 
r ~ 0.4_Ep. The perturbative physics of the adiabatic regime 
is approached only as F exceeds E p , leaving an intermediate 
crossover region for 0AE p < T < E p where the polaron gets 
progressively undressed. The phonon configuration strongly 
deviates in the crossover region from a simple superposition 
of just one or two coherent states. 



traditional antiadiabatic regimes is therefore quantita- 
tive rather than qualitative, both being described by the 
same fermionic interacting resonant-level model at ener- 
gies below ujq. For this reason, we have lumped the two 
regimes into one in Fig. [T] In terms of the bare model 
parameters, the extended antiadiabatic regime persists 
up to r ~ 0.4£" p , breaking down as T cS approaches uiq in 
magnitude. The extended antiadiabatic regime and the 
adiabatic one are separated in turn by a true crossover 
region for 0AE p < T < E p , where the polaron gets pro- 
gressively undressed. The phonon configuration strongly 
deviates in this region from a simple superposition of just 
one or two coherent states, as we show by explicit calcu- 
lations. 



B. Plan of the paper 

After introducing the Hamiltonian and its symmetries 
in Sec. [TTl we proceed in Sec. IIIII to a preliminary dis- 
cussion of certain limits where analytic insight can be 
gained. These include the traditional adiabatic and an- 
tiadiabatic limits, the extended antiadiabatic regime pro- 
posed in this paper, and the limit of large detuning. Sec- 
tion llVl presents in turn a systematic study of all coupling 
regimes using Wilson's NRG. To this end, we begin with 




FIG. 2: (Color online) The physical system under considera- 
tion. A localized level with energy e d is tunnel coupled with 
amplitude t to a band of spinless electrons, and is simultane- 
ously coupled by a dimensionless displacement coupling A to 
a single vibrational mode of frequency uuo. Depending on the 
valence of the level, the vibrational mode experiences a shifted 
harmonic potential with a relative displacement of Ab = A. 



a brief introduction of the method in Sec. IIV A[ followed 
by an extensive investigation of the key quantities of in- 
terest: the renormalized tunneling rate T e g, the map- 
ping onto an effective low-energy interacting resonant- 
level model in the antiadiabatic and the extended antia- 
diabatic regimes, charging of the level, and the phonon 
distribution function defined in Eq. (|35[) . We conclude in 
Sec. I VII with a summary of our results. Some technical 
details are deferred to two appendices. 



II. THE MODEL AND ITS SYMMETRIES 

The system under consideration is shown schematically 
in Fig. [5] It consists of a single localized electronic level 
$ with energy e d , tunnel coupled to a continuous band 
of noninteracting spinless electrons which we denote by 
c\. The level is simultaneously coupled by displacement 
to a localized vibrational mode (phonon) , as modeled by 
the Hamiltonian 

n = j2^ c k c k + ^J2{ dic k +R - c -} + e d^ d 

k V i V k 

+Lo tfb + \uj Q (dU-N )(b i < + b). (1) 

Here, W creates a local Einstein phonon that oscillates 
with frequency u; , t is the tunneling matrix element 
between the level and the Wannier state closest to the 
molecule, and N is the number of lattice sites. The di- 
mensionless coupling A measures the relative displace- 
ment of the vibrational mode between the configurations 
where the level is empty and occupied. It serves as a 
faithful measure for the strength of the electron-phonon 
coupling, with A C 1 (1 < A) corresponding to weak 
(strong) interactions. The parameter No can be thought 
of as fixing the reference charge of the level. It can be 
formally eliminated by shifting the bosonic mode accord- 
ing to B = b — XNq, which has the effect of renormalizing 
the level energy from its bare value e d to e d with 



(2) 



The conversion from b to B also generates the constant 
term — iV A 2 o;o, which uniformly shifts the spectrum of 
the Hamiltonian. Although the inclusion of Nq adds no 
richness to the thermodynamics of the model, it provides 
a useful tuning parameter for exploring the low-energy 
state of the system, as will be demonstrated later on. 

In the following we shall consider a particle-hole sym- 
metric band, namely, it is assumed that the wave num- 
bers can be grouped into distinct pairs k and k', such 
that efc' = — £fc for each pair of momenta. It is easy to 
verify that the combined transformation 



d) -> d, ci -> - 



Cfc' 



and b -> -b - A(l - 2N ) (3) 



leaves the Hamiltonian of Eq. (JXJ) unchanged under these 
terms, apart from the substitution^ 2 - 



l d = 2A 2 w (l - 2AT ) 



(4) 



It therefore suffices to study the domain e d > with 

e* d = A 2 c (l - 27V ), (5) 

while the complementary domain e d < e* d is accessible 
via the particle- hole transformation of Eq. ([3]). In par- 
ticular, occupancy of the level rid(e d ) = (d^d) obeys the 
symmetry relation 



^(e3-Ae d ) = l-n d (e*+Ae d ), 



(G) 



independent of all other model parameters, the temper- 
ature included. 

Other than the conduction-electron bandwidth D, the 
Hamiltonian of Eq. ([T]) features four basic energy scales. 
These include the bare vibrational frequency u>o, the po- 
laronic shift E p = A 2 wo, the detuning energy Ae d = 
€ d ~ e d> anc ^ the hybridization width T = irpot 2 . Here 
Po is the conduction-electron density of states at the 
Fermi levels Of particular interest is the point Ae d = 0, 
when the (renormalized) level lies at resonance with the 
Fermi energy. The relevant low-energy scale in the prob- 
lem is conveniently defined 2 ^ in this case from the zero- 
temperature charge susceptibility evaluated at e* d : 



efit 



1 



with 



dn d 
de. 



(8) 



For A = 0, the low-energy scale T e g so defined coin- 
cides with the bare hybridization width T. Various ther- 
modynamic properties associated with the level, e.g., its 
occupancy and its contribution to the electronic specific 
heat, reduce in the wide-band limit to exclusive functions 
of e d /Y and T/Y, where T is the temperature. A nonzero 
A modifies this picture both qualitatively and quantita- 
tively. For example, occupancy of the level at T = 



(a) <b> 



(b) G = 



(c) 



(d) 



(e) 



- A.Q)„N h 



o 



0(X 3 ) 



FIG. 3: Diagrammatic representation of (a) the phononic ex- 
pectation value (b) and (b) the phononic Green's function 
G(iv„). Here, single (double) wiggly lines represent the bare 
(dressed) connected phonon propagator G(iv n ), whereas sin- 
gle (double) lines with an arrow denote the bare (dressed) 
d-electron propagator (bare/dressed with respect to A; both 
electronic propagators are dressed with respect to t alone). 
The connected phonon propagator G(iv n ) obeys a stan- 
dard Dyson equation (d), with a matrix self-energy E whose 
second-order expansion in A is specified in panel (e). Interac- 
tion vertices are marked by a cross. 



is no longer given for large A by an exclusive function of 
Ae d /r e flf, but rather depends, as we shall show, on vastly 
different energy scales. Our goal is to conduct a system- 
atic study of all coupling regimes for nonzero A, focus- 
ing primarily on large A. To this end, we shall combine 
analytical considerations with Wilson's renormalization- 
group (NRG) method^^ 



III. LIMITING CASES 

We begin with a preliminary discussion of certain lim- 
its where analytic insight can be gained. These include 
the traditional adiabatic and antiadiabatic limits, the 
proposed extended antiadiabatic regime where luq < Y 



(7) but r cff -c wo, and the limit of large detuning, Y <C |e<j| 



A. Adiabatic limit 

Commencing with the limit of small A, we apply ordi- 
nary perturbation theory in the electron-phonon coupling 
A, postponing for the moment the question of its range of 
validity. The basic quantity of interest is the Matsubara 
phonon propagator 



G{iv n 



G b tf{w n ) G bb {iv n ) 
G b f b f(ip n ) G b t b (iv n ) 



(9) 
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where v n — 2irn/f3 are the bosonic Matsubara frequen- 
cies, /3 — l/k^T is the reciprocal temperature, and 



G AB {w n ) = - / (T T A(T)B(0))e l ^ T dT. 



(10) 



Using the diagrams specified in Fig. [3j one obtains the 
formally exact relations 



(b) = -X [(dU) - N ] 



(11) 



and 



G aj (iv n ) = G ai (iv n ) ~ 5 nfi p\ 2 [(dU) ~N ]\ (12) 

where G{iv n ) is the connected phonon propagator, de- 
fined as the sum of all connected phonon diagrams. The 
connected propagator has the conventional representa- 
tion G" 1 = [G^]- 1 ~ S, where 



iv n - loo 











-iv n - w 



(13) 



is the unperturbed phonon Green's function and 2 is the 
self-energy matrix. Since the electron-phonon interaction 
in Eq. ([T]) involves only the combination b + b\ the self- 
energy matrix takes the general form 



with 7 = X 2 uj 2 o/{itT 2 ) and A = \ 2 ul/{irT). The poles of 
G(z) can now be identified with the zeros of det{G _1 (z)}, 
which yields 



,w / Aw \ , L 

z± = - l -\—) ± "°v 1 



Aw 

r 



tt r 



(20) 

Several conclusions can be drawn from Eq. (f20|) . First, 
there are two distinct parameters that control the per- 
turbative expansion: E p /T and Xwq/T. Depending on 
the magnitude of A either parameter can be the largest, 
with Ep/T dominating for 1 <C A. Second, both parame- 
ters must be small in order for a weakly damped phonon 
mode to persist. A crossover to an overdamped phonon 
occurs as soon as either of these parameters becomes of 
order unity, signaling a qualitative change in the under- 
lying physics and the breakdown of perturbation theory. 
Third, coupling to the electronic level softens the phonon 
frequency according to 



ejeff _ 1 _J_f^o\ 2 E p 

ujo y 7r 2 \ r y 7T r 



(21) 



For 1 C A we therefore conclude that the perturbative 
physics of the adiabatic limit breaks down as soon as the 
polaronic shift exceeds T. 



S(ii/„) = a(iu n ) 



1 1 
1 1 



(14) 



which depends on a single scalar function a(iv n ). Settling 
with second order in A and analytically continuing to real 
frequencies, a(uj + irj) is given for T — > in the wide-band 
limit by 



a (to +irj) 
where 



2in 



^[/_( w )-J*(-o;)-J + H+4(- W )], 

(15) 



1+ («) = - [ln(-e d + tT) - In(w -e d + iT)] (16) 

U! 



and 
J_( W ) 



1 



ui + 2iT 



[iir + h\(e d + iT) - ln(w - e d + iT)} 



(17) 

Focusing on the resonance condition e d = and expand- 
ing in powers of the frequency— a(uj + iff) reads as 



a{ui + irj) 
resulting in 



\ 2 / i 2 \ 2 , i 2 



(18) 



(1 + 47)0; — uj + A 
i^uj + A 



ijuj + A 

(ij—l)uj — ujo + A 

(19) 



B. Antiadiabatic limit 

When r is sufficiently small as compared to ljq, all 
phonon excitations are frozen out as the temperature 
is decreased below uiq. Only fermionic excitations re- 
main active at such low energies, reflecting the fact 
that the relevant electronic motion is far slower than 
the phonon vibrations. To derive the effective low- 
energy Hamiltonian for Ae d — 0, it is useful to first ap- 
ply the Lang-Firsov transformational H' = U^TiU with 
U = e -A(ft t +ft)(d t d-Ar ) j which conv e r ts Eq. © into 



-UJ 



k 

(22) 

Here we have set Ae d = and omitted the global energy 
shift AE = -X 2 ujoN§. The effect of the Lang-Firsov 
transformation is to eliminate the displacement interac- 
tion term at the expense of attaching the exponentials 
e ±A(6+-fc) £ Q £ ne tunneling amplitude t. The transformed 
Hamiltonian is invariant under the particle-hole transfor- 
mation 



-> -c 



k'> 



and b — > —b, 



(23) 



which comes in place of Eq. . 

Provided T is small enough, the (transformed) phonon 
is frozen in its unperturbed ground state configuration 
&t& = at energies well below ojq. The effective Hamilto- 
nian at such low energies is purely fermionic, and may 
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generally contain all possible local Hamiltonian terms 
that are invariant under the particle-hole transformation 
<$ — s- d, c]. — > —c k ,. A rigorous derivation of the effec- 
tive low-energy Hamiltonian requires a systematic inte- 
gration of all high-energy excitations, including a succes- 
sive elimination of all discrete phonon excitations. This 
turns out to be a difficult task. Nevertheless, for T <C wo 
one can settle with a single-step elimination of all excita- 
tion energies exceeding ujq using a Schrieffer- Wolff-type 
transformation.-- Deferring the details of the derivation 
to Appendix IA1 we quote here only the end result. 

To second order in t, one is left with an effective inter- 
acting resonant-level model 27,28 (IRLM) of the form 

restric , rcstric 

n cS = J2 + -7= E l rftc * + H - c -} ( 24 ) 

k v k 



Here, 



N 

c l c k> : = 



(dU-1/2) J2 :c K<: 



c i c k> 



9(~£k)5k.k' stands for normal or- 



dering with respect to the filled Fermi sea, while the sym- 
bol XTfc° 8 " C comes to indicate that the summation over k 
is restricted to momenta such that |ejt| < ujq (N' < N be- 
ing the number of such k points). The coupling constants 
entering Eq. (|24[) are given by 



t c ff = te 



-A 2 /2 



(25) 



and 



9t 2 

U eS = —e- x [Ei(A 2 ) - 2 ln(A) - 7] , (26) 

UJQ 

where Ei(x) is the exponential integral function 2 — and 
7 « 0.5772 is Euler's constant. In the limits where A is 
either small or large as compared to one, U e g takes the 
simplified forms 



U, 



eff 



2i 2 A 2 /wo, A<1, 
2t 2 /(A 2 w ), 1 « A, 



(27) 



as follows from the corresponding asymptotes of the ex- 
ponential integral function. 

One can borrow at this point known results for 
the IRLM in order to extract the renormalized level 
width r o ff. Specifically, it is known from perturbative 
renormalization-group calculations^ that 



r cff = d. 




l/(l+2p C/ o ff) 



(28) 



where _D c ff is the effective bandwidth and T = npot^g. 
Inserting Eq. (f25j) for t c e and setting D c ff ~ ujq yields 



T c ff = w ( — e 



—A* 



l/(l+2p C/ oft ) 



(29) 



It should be stressed that Eqs. (|25|). (|26j). and (J29J) are 
restricted to T <C M)> when the exponent that appears 
in Eq. (|29|) hardly deviates from one. As we show be- 
low, the reduction to an effective IRLM is far more gen- 
eral, though, and applies to all parameter regimes where 
L e ff < M). 



C. Extended antiadiabatic regime 

Once F exceeds uiq, one looses the hierarchy of en- 
ergy scales underlying the derivation of the effective low- 
energy Hamiltonian of Eq. ([Ml) . Nevertheless, the gen- 
eral form of the low-energy Hamiltonian can still be de- 
duced for T c ff <C w on physical grounds. Whenever 
F c ff <C Wo, the renormalized electronic motion remains 
sufficiently slow for the phonon to efficiently respond to 
the individual tunneling events. Hence, one can expect a 
rather well-defined bosonic mode to persist, with a mod- 
ified frequency that remains close in magnitude to ujq. 
As in the strict antiadiabatic limit, the bosonic mode is 
frozen in its ground-state configuration as the tempera- 
ture is lowered below u>o, leaving only purely fermionic 
excitations active at such low energies. With increasing 
F, the microscopic details of the bosonic mode and its 
associated polaron will progressively deviate from their 
F — > forms, yet the effective low-energy Hamiltonian 
remains purely fermionic as long as T e ff <wo. As we ar- 
gue below, the general form of the resulting Hamiltonian 
is essentially dictated by symmetry considerations. 

To illustrate this point, let us begin with the trans- 
formed Hamiltonian of Eq. (|2"2"|) . As indicated above, the 
effective fermionic Hamiltonian may generally contain all 
possible local Hamiltonian terms that are invariant under 
the particle-hole transformation e?t — > d, c\ — > — c k ,. Of 
all possible terms in this category, only the local tunnel- 
ing term t is relevant in the renormalization-group sense, 
while the local contact interaction 



u(dU- 1/2) -oW 



(30) 



k,k' 



is marginal. All other local terms permitted by sym- 
metry are formally irrelevant, leaving us with the ef- 
fective IRLM of Eq. (|24l) . A nonzero Ae d relaxes the 
requirement of particle-hole symmetry, which augments 
the Hamiltonian of Eq. (IM1) with two additional terms: 



£cff 



dU- 



VeS 

N 



E 

k,k' 



'■ c k C k' 



(31) 



Note that although the form of Eq. (|24l) is dictated 
by symmetry considerations, these considerations alone 
do not suffice to fix the values that the couplings t e g and 
U e s acquire. These couplings must generally be extracted 
from the low-energy spectrum of the original electron- 
phonon Hamiltonian, as will be done later on using the 
NRG. As we shall show, the condition T c ff <C uj en- 
compasses for I <§; A all values of F up to F ~ QAE p . 
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Furthermore, U e s remains well described by Eq. (|26|) for 
most of this range, whereas t e g rapidly exceeds Eq. (j2"5)l 
as soon as T approaches loq. 



D. Large detuning 

One limit where perturbation theory in t is guaranteed 
to apply is that of large detuning, Ae d S> T (equivalently 
— Ae d T). In this case Ae d is sufficiently large to 
assure that the electronic level remains nearly empty at 
zero temperature, providing a suitable starting point for 
a perturbative expansion in t. In the following we focus 
on the zero-temperature occupancy of the level, rid(e d ), 
which is conveniently computed from the derivative of 
the ground-state energy with respect to e d . 

To obtain the correction to the ground-state energy it 
is useful to start from the transformed Hamiltonian of 
Eq. (j2"2")l , which is augmented for Ae^ ^ by the Hamil- 
tonian term Ae d <fid. For t = and Ae d > 0, the ground 
state of H' is given by the product state of the filled Fermi 
sea with an empty level and the empty phonon state. To 
second order in t the ground-state energy acquires the 
correction 



distribution function 



SE gs = — e 

7T 



r _ x2 ~ a 2 " 



n=0 



In 



Ae d + D + nuj 



(32) 



where we have assumed a symmetric rectangular density 
of states for the conduction electrons: pie) — pod(D— |e|). 
Straightforward differentiation of Eq. (|3"2"|) with respect to 
ed yields then the level occupancy 



7r £ n! V At, 4 



n=0 



d + nuj 
f 



(33) 



D + Ae , + nu>o 



which properly reduces for A — > to the noninteracting 
result 



n d (e d > 0) = 



I 

AT" 



1 



D 



(34) 



In the wide-band limit, D — > oo, the second term drops 
out in the parentheses of Eq. (|33p . 

Note that Eq. (|33l) was derived under the strict condi- 
tion that Ae d 3> T. Below we present NRG results that 
suggest a broader range of validity of Eq. (|3"3")l . down to 
Ae d ~ r off . Since r off < T for 1 < A and r <C E p , this 
implies a far greater range of validity of perturbation the- 
ory in t. 



E. Phonon distribution function 

Our discussion thus far was restricted to electronic 
properties. Another quantity of interest is the phonon 



P(n) = (\n)(n\), 



(35) 



which contains direct information on the state of the 
phonon. The phonon distribution function has distinct 
characteristic forms in the extreme adiabatic and antia- 
diabatic limits, which we next derive. The transition be- 
tween these two limiting forms indicates the undressing 
of the polaron upon going from the antiadiabatic to the 
adiabatic regime. 

In the perturbative adiabatic regime, the phonon is 
too slow to respond to the successive tunneling events, 
hence it samples only the time-averaged occupancy of 
the level. From the standpoint of the phonon it therefore 
experiences the effective Hamiltonian 



■Hphonon = ojob^b + XuJa{n d - N ) (6 f + b) 



(36) 



describing the average displacement A — > A = \(n d — N a ). 
At T — the phonon is thus frozen in the coherent state 
£> | A) = — A | A), resulting in 



adiabatic 



(n) 



A 



2» 



(37) 



For the particular case where Ae d = and Nq = 0, one 
has that A = A/2. 

In the antiadiabatic limit, it is advantageous to con- 
sider first the ground state \tp' gs ) — C^lV'gs) °f the 
transformed Hamiltonian H' of Eq. ([22]) . Here U — 
e -A(f) t +6)(d t d-A'o) j s -j-^g canon i ca j transformation relat- 
ing H and H' . As discussed in Sec. IIII B[ for r <C ojq 
the transformed phonon is effectively frozen at T = in 
its unperturbed ground state Wb — 0. This means that 
\iptj) is well approximated by a product state of the form 



Ks) = 10) ci ® |n = 0) ph , 



(38) 



where \4>)ei pertains to the electronic degrees of freedom. 
The ground state of the original Hamiltonian l^gs) can 
now be obtained by applying the transformation U to 
iV'gg), resulting in 

\ip' ss ) = |0o)d ® |A ) P h + |0i)oi ® |Ai) ph . (39) 

Here |^o)ei and \4>i} e i, respectively, are the projections 
of the electronic state \4>) e \ onto the n d — and h d — 
1 subspaces, while |Ao) p h and |Ai) p h are the phononic 
coherent states with Ao = AA^o and Ai = X(Nq — 1). 
Accordingly, the phonon distribution function assumes 
the form 



A 

Panti-adiabatic(w) = ( 1 - Tl d ) -^7~e~ X ° + n d -^T 



2;i 



2n 



III 



In particular, for Nq = Ae d = 0, 



P, 



ant i— adiabatic 



(n) 



(40) 



(41) 
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IV. SYSTEMATIC STUDY OF ALL COUPLING 
REGIMES 

A. The Numerical renormalization group 

To treat the Hamiltonian of Eq. (JXJ) for arbitrary 
coupling strengths, we resort to Wilson's numerical 
renormalization-group metho d 19 ' 20 (NRG). The NRG is 
a powerful tool for accurately calculating equilibrium 
properties of arbitrarily complex quantum impurities. 
Originally devised for treating the single-channel Kondo 
Hamiltonian^ this nonperturbative approach was suc- 
cessfully applied over the years to numerous impurity 
models and setups^ At the heart of the approach is a 
logarithmic energy discretization of the conduction band 
about the Fermi energy, controlled by the discretization 
parameter A > 1. Using an appropriate unitary transfor- 
mation^ the conduction band is mapped onto a semi- 
infinite chain with the impurity coupled to its open end. 
The Nth link along the chain represents an exponentially 
decreasing energy scale Da ~ A _Ar / 2 , with the contin- 
uum limit recovered for A — > 1 + . The full Hamiltonian 
of Eq. ([T]) is thus recast as a double limit of a sequence 
of dimensionless NRG Hamiltonians: 



H = 



lim Km {d a A- (N - 1)/2 H n \ 



(42) 



with D A = D(l + A- 1 )/2 and 



N 



A— 



+Aw (Sd - N ) (fet + b) 



N-l 



x;a-^„{/+ +1 /„+h.c} 



(43) 



Here, id = e^/ 'Da and ljq = ujq/Da are the dimensionless 
energy level and vibrational frequency, respectively, while 
t is related to the tunneling matrix element through 



The coefficient 



t 



A^ = 



A + l 

2(A-1) 



In A 



(44) 



(45) 



is required to account for the energy discretization used 
in the NRG,— and can be viewed as accelerating the 
convergence to the A — » 1 + limit. The prefactor A^ -1 ^ 2 
that appears in Eq. (|4"3"|) comes to ensure that the low- 
lying excitations oiHA are of order one for all N. 

Physically, the shell operator /J represents the local 
conduction-electron state to which the level is directly 
coupled by tunneling. The subsequent shell operators /t 
correspond to wave packets whose spatial extent about 
the level grows roughly as A"/ 2 . Details of the band 
are encoded in the hopping coefficients £„, obtained from 



suitable integrals of the density of states^ 2 - Throughout 
the paper we assume a relativistic dispersion relation, 
corresponding to the symmetric rectangular density of 
states p(e) = p a 9(D~\e\) with p a = 1/(2D). This simpli- 
fied form of p(e) affords an explicit analytical expression 
for £ n r^ which rapidly approaches one with increasing n. 

A key ingredient of the NRG is the separation of energy 
scales along the Wilson chain, which enables an iterative 
diagonalization of the sequence of finite-size Hamiltoni- 
ans T~Ln- Starting from a core cluster that consists of the 
local degrees of freedom d\ b\ and /g, the Wilson chain 
is successively enlarged by adding one site at a time. Di- 
agonalization of Hn+i proceeds from the knowledge of 
the spectrum of Hat by means of the NRG transforma- 
tion 



n N+1 = Vxn N + 6v {fl+Jn + h.c.} 



(46) 



In this manner, one can track the evolution of the finite- 
size spectrum as a function of N. The approach to a fixed 
point is signaled by a limit cycle of the NRG transforma- 
tion, with "Hat+2 and Hat sharing the same low-energy 
spectrum. 

The above procedure could, in principle, be applied 
to any set of hopping matrix elements along the chain. 
However, practical considerations prove far more restric- 
tive, as it is numerically impossible to keep track of the 
exponential growth of the Hilbert space with increasing 
N. In practice only a limited number of states can be 
retained at the conclusion of each NRG iteration, which 
is where the separation of scales along the chain comes 
into play. Due to the exponential decrease of the hop- 
ping terms, one can settle with retaining only the lowest 
N B eigenstates of Ha when constructing the low-energy 
spectrum of Hn+i- The NRG eigenstates so obtained 
are expected to faithfully describe the spectrum of Hn 
on a scale of Das = -DaA - ^ -1 )/ 2 , corresponding to the 
temperature Taj ~ Da- Thus, three distinct approxima- 
tions are involved in the NRG algorithm when applied 
to the Hamiltonian of Eq. (JTJ: (i) Discretization of the 
conduction band, controlled by the parameter A > 1; (ii) 
A finite-size representation of the bare bosonic spectrum, 
controlled by the number A^ of bare bosonic states kept 
(we use the states where wb = 0, • • • , N\, — 1); (iii) Trun- 
cation of the Hilbert space at the conclusion of each NRG 
iteration, controlled by the number N s of states retained. 
Each of these three approximations can be systematically 
improved by varying A, N^, and N s . All data points pre- 
sented in this paper were obtained for A = 2, while the 
number of states retained were either N\, = 1600 and 
A^ = 4000 or N h = 3000 and N s = 8000. Explicit values 
are quoted in the relevant figure captions. 

We emphasize that the total number of electrons, or Q 
in the notation of Ref. l3ll . is the only conserved quantity 
one can exploit in the iterative diagonalization of the 
sequence of NRG Hamiltonians for our problem. Since 
each additional site can either be empty or occupied, it is 
straightforward to keep track of the associated quantum 
number using the algorithm detailed, e.g., in Ref. l3lL 
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FIG. 4: (Color online) The ratio F eff /r vs R = r/_B p , for 
ojo/D — 0.0025 and different strengths of the electron-phonon 
coupling A. Here, A takes the values 3 (red), 4 (green), 5 
(blue), 6 (orange), and 8 (maroon). Symbols depict the NRG 
data. Solid lines show the ratio r e mp/r, obtained using the 
empirical formula r cmp = r exp[— A 2 J-(R)] with J-(R) given 
by the solid line in Fig. [5] The gray shaded area covers the 
regime where the condition T c ff = ojq is met for all values 
of A displayed. Inset: The ratio r e ff/r em p for all data points 
displayed in the main panel. The empirical formula and NRG 
data agree to within a factor of 2.3 for all values of A and all 
ratios T/E p displayed (the dotted line marks the value 0.44). 
NRG parameters: A = 2, iV s = 8000, and N h = 3000. 



B. Renormalized tunneling rate r e fi 

We begin our discussion with the renormalized tunnel- 
ing rate r e ff, denned by Eqs. (0 and (|SJ) with \c eval- 
uated for T —> 0. Figure 0] shows the dependence of 
r e g on r for ujo/D = 0.0025 and different strengths of 
the electron-phonon coupling A. As expected, the ratio 
r e ff/r varies by orders of magnitude upon going from 
small to large values of T when A is large. For E p <C T, 
one essentially recovers the bare tunneling rate T, while 
for r <C wo there is an exponential suppression of the 
tunneling rate according to r e gf/r = e _A . Remarkably, 
the crossover between these two limits follows an approx- 
imate scaling form, as demonstrated in Fig. [5] Plotting 
A~ 2 ln(r/r c ff ) as a function of R = T/E p for the different 
values of A, all data points approximately collapse onto 
a single curve. The collapse is particularly good in the 
range R < 0.2, and gradually degrades for larger values 
of R. Note that E p /D equals 0.16 for A = 8, hence some 
values of T become of order the bandwidth for R ~ 10. 

The quality of the approximate scaling form (which, 
as we show below, is not an exact scaling function) can 
be appreciated by extracting an empirical function J-(R) 
such that ^(R) ~ A~ 2 ln(r/r e ff), and comparing the cal- 
culated values of r c ff to the empirical formula 

r cmp = rexp[-A 2 ^(i?)]. (47) 

The empirical formula, depicted by the solid lines in 
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FIG. 5: (Color online) The function A~ 2 ln(r/T cff ) vs R = 
r/E p , for ujq/D — 0.0025 and different strengths of the 
electron-phonon coupling A. The solid line is an empirical 
curve that defines the function J-(R), which is used in Fig. [4] 
to compute F cmp — A GXp {-\ 2 T(R)]. All NRG parameters 
are the same as in Fig. [4] 

Fig. 21 well agrees with the calculated values of r e ff over 
many orders of magnitude. Deviations are confined to 
within a factor of 2.3 (see inset of Fig. 2]), which is quite 
remarkable considering the enormous variation in r e ff as 
a function of both A and T. As for the function T(R), we 
define it by the solid line in Fig. [5] Obviously, there is 
some arbitrariness in the way ^(R) is fixed, particularly 
for R > 0.2 where scaling degrades. Indeed, it is in this 
parameter regime that the deviations between r emp and 
r o ff are typically the largest. Still, Eq. (T4"T|) provides a 
useful formula for the renormalized tunneling rate, suc- 
cessfully interpolating between the extreme adiabatic and 
antiadiabatic limits. Finally, we note that the condition 
r c ff = cj is met for R in the range 0.53 < R < 0.59 (gray 
shaded area in Fig. 2]) for all values of A displayed, hence 
the two scales remain well separated up to R ~ 0.4. 



C. Extended antiadiabatic limit: Mapping onto the 

IRLM 

Next we focus on the resonance condition Ae^ = and 
examine in greater detail the regime where T c g <C lo^. 
For the case of interest where 1 <S A, this condition cor- 
responds to r <C E p , which encompasses both the antia- 
diabatic limit r < and an extended region where T 
exceeds ujq and yet r e ff <c ujq- As discussed in Sec. IIII Cl 
we anticipate for such couplings that the system is de- 
scribed at energies below ojo by an effective IRLM with 
the tunneling amplitude t c ff and the local Coulomb repul- 
sion U c s- Using the NRG level flow, we have confirmed 
this physical picture. In the intermediate energy regime 
r e g <C Dn < wo, the finite-size spectra consistently re- 
duced to that of a weakly coupled IRLM, proving the 
validity of the extended antiadiabatic regime. The cou- 
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FIG. 6: (Color online) The tunneling amplitude t e ff that 
enters the effective IRLM at energy D e ff/u;o = 0.146, plot- 
ted vs R = T/E p . Here ujq/D = 0.0025, while D e s equals 
Dn=22 with A = 2. The values of A are 3 (red), 4 (green), 5 
(blue), and 6 (orange). The remaining NRG parameters are 
N s = 8000 and Nb = 3000. For clarity, i e ff was normalized by 

— A 2 II 

^analytic = te , which is the effective tunneling amplitude 
obtained for r<wo- Inset: A scaling plot of t e ff/tanalytic vs 
T/loq. The solid line shows the parabola f(x) — l + 0.9x + x 2 . 



pling constants that enter the effective Hamiltonian can 
be read off from the NRG spectra using the procedure 
outlined in Appendix [Bj Our results are summarized in 
Figs. © and 

Figure [6] displays the effective tunneling amplitude t e ff 
at energy D c s/ll>q — 0.146 for several values of A = 2, 4, 
5, and 6. For clarity, all curves have been normalized by 
^analytic = t e ~ A which accounts for the main Gaussian 
dependence of t e ff on A. As can be seen in the inset, all 
data points for t e ff /tanalytic collapse onto a single curve 
when plotted versus T/luq, at least for values of T up to 
a few times loq. Thus, the effective tunneling amplitude 
acquires the empirical scaling form 



t e g = t ■ exp 



A 



2 1 



/(T/wo), 



(48) 



where f(x) is well fitted by the parabola f(x) ~ 1+0. 9x+ 
x 2 (depicted by the full line in the inset of Fig. [6j . Three 
points are noteworthy. First, t c ff grows quite rapidly with 
r/wrj, increasing by a factor of 5 in the limited range 
covered by the inset of Fig. [6] Second, given the apparent 
scaling of t e g /t a naiytic with r/wrj, it is evident that the 
ratio r c ff/r cannot be an exclusive function of F/E p as 
suggested by the scaling plot of Fig.[S] Lastly, we cannot 
reliably extract t e g for larger values of T/ojq since t e ff and 
t/eff no longer represent well-separated energy scales (for 
technical details, see Appendix |B|) . 

The effective Coulomb repulsion U c g at energy 
D e s/uJo = 0.146 is shown in turn in Fig. after divi- 
sion by its asymptotic weak-tunneling form t/ ana i yt i C — 
2R/(ttp ) [see Eq. d^TJ) with 1 < A]. Note that in con- 
trast to t ff, which becomes numerically inaccessible for 
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FIG. 7: (Color online) The local Coulomb repulsion U e s that 
enters the effective IRLM at energy D e g/cjo = 0.146, plotted 
vs R — F/Ep. The values of A are 5 (red), 10 (green), and 20 
(blue). All other parameters are the same as in Fig. [5] For 
clarity, U e a was normalized by ^analytic = 2R/(ivpo), which is 
the effective Coulomb repulsion for r ujo. Inset: A scaling 

plot Of C/eff/t/analytic VerSUS R/\f\. 



8 < A, the effective Coulomb repulsion U e ff can be ac- 
curately computed for values of A well above 10. Sur- 
prisingly, the weak-tunneling expression that was derived 
strictly speaking for r <C ujq remains quite accurate (to 
within 10%) even for T/ujq as large as 100 when A = 20. 
Hence, the main source of T dependence stems from 

[^analytic which Scales as T/E p . Similar to £ e ffAanalytic 

also t/eff /[^analytic appears to follow an approximate scal- 
ing form, this time with the scaling variable r/(A 2 5 wo) 
(see inset of Fig. [7]). However, the quality of the data 
collapse and the variation in t/ e ff/t/ a naiytic are far more 
restricted than for t ff/t a nai y tic- 

From the discussion above it is clear that the regimes 
r < wo and ujq < T <C E p share the same qualitative 
physics, both being described by the same IRLM at en- 
ergies below bjQ. The distinction between the two regimes 
is mainly quantitative, as encoded in the effective model 
parameters i e ff and U e g. The rather rapid departure of 



t G ff from its asymptotic weak-tunneling form t e 



-A72 



re- 



flects its extreme sensitivity to even small deformations 
of the polaronic mode. 



D. Charging of the level 

Up until now, our discussion was restricted to Ae^ = 0. 
Next we consider nonzero detuning and examine the 
charging properties of the d level. Besides being of inter- 
est on its own right, the charge of the level is intimately 
related at T — to the conductance of the molecular 
bridge depicted schematically in Fig. [TTJ We address the 
latter setup in detail in Sec. IVl 

Figure [8] shows the level occupancy rid versus Ae^ > 0, 
for T 4 and A = 6. The level occupancy in the com- 
plementary regime Ae d < is obtained from the sym- 
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FIG. 8: (Color online) The level occupancy n d vs Ae^ = 
e d - e* d > 0, for uj /D = 0.0025, T -> 0, A = 6, and different 
values of J? = T/E p (indicated by arrows). The level occu- 
pancy in the complementary regime Ae d < is obtained from 
the symmetry relation of Eq. |(5J. With increasing F > E p , 
charging of the level approaches the conventional noninter- 
acting curve, depicted by the dotted line. With decreasing 
r < E p , the noninteracting shape rapidly deforms into a 
double-step structure governed by distinct mechanisms for 
Ae d < r c fr and F c ff <C Ae d . For Ae d < r e ff, one recovers 
a strongly renorrnalized noninteracting form with V — > r o ff. 
For r e ff <C Ae d , charging is well described by a simple per- 
turbative expansion in t. This is demonstrated in the inset, 
where each of the occupancies of the main panel is compared 
to the perturbative expression of Eq. (|33[1 (dashed line). NRG 
parameters: A = 2, iV s = 4000, and iV b = 1600. 



0.1 


0.1 


0.1 


0.1 

0, 



.0.49 



• 0.36 



R = 0.1 



R = 0.3 











irrni 



R = 0.6 



R = 2.0 



* ■ ■ 







10 15 20 25 
n 



FIG. 9: (Color online) The phonon distribution function 
P(n) = (\n)(n\), for u /D = 0.0025, X = 4, N = 0, T -> 0, 
and different values of R = T/E p . The detuning energy Ae d 
is set to zero such that the electronic level is half filled. For 
R <C 1, exemplified by R = 0.1 in the upper panel, the distri- 
bution function P(n) approaches the double-peak structure 
of Eq. (|4ip . depicted by the dotted line. In the opposite 
limit where R is large, exemplified by R = 2 in the lower 
panel, the distribution function P(n) approaches the single- 
peak structure of Eq. (|37p . displayed by the dotted line. NRG 
parameters: A = 2, iV s = 4000, and N h = 1600. 



which is insufficient for a fully developed shoulder to be 
seen. 



metry relation of Eq. ([5]). With increasing r > E p , 
charging of the level approaches the conventional non- 
interacting curve, depicted by the dotted line. With de- 
creasing r < E p , the noninteracting shape rapidly de- 
forms into a charging curve governed by two distinct 
mechanisms: (i) a strongly renorrnalized noninteracting 
form with T — > T c tf, applicable up to Ae d ~ r c ff, and 
(ii) a simple perturbative expansion in t, applicable for 
T c ff <C Ae d . The latter mechanism is demonstrated in 
the inset, where each of the curves of the main panel 
is compared to the perturbative expression of Eq. ([55)1 
(dashed line). 

The combination of these two mechanisms gives rise to 
a distinctive shoulder in rid, which interpolates between 
the limits where Ae^ ~ T c ff and Ae d ~ T. The height 
of the shoulder decreases with decreasing T, approaching 
T/(ttE p ) when r <C E p . This result can be understood 
from the fact that the summation over n in Eq. (|33[) 
samples mainly the regime where n ~ A 2 when A is large, 
hence the denominator Ae^ + uujq is effectively replaced 
with Atd + A 2 ojo ~ E p . This argumentation breaks down 
as Aed approaches T e g, when higher order terms become 
exceedingly more important. Another effect of decreasing 
r is the opening of an exponential separation between T 
and T e ff, which sets the lateral extent of the shoulder 
when plotted versus Ae d /T e g. For example, T and T e g 
are separated by one order of magnitude for R = 0.6, 



E. Phononic distribution 

To gain direct information on the state of the phonon, 
we next consider the phonon distribution function P(n), 
defined in Eq. pg]). Figure [9] depicts P(n) for T ->• 
using A = 4, 7V = 0, and different values of R = T/E p . 
The detuning energy Ae d (which itself depends on A*o 
through e* d ) is set to zero, corresponding to resonance 
condition. For small values of T, exemplified by R = 0.1 
in the upper panel, the distribution function P(n) ap- 
proaches the double-peak structure of Eq. (|4"T]) , which is 
plotted for comparison by the dotted line. Although T 
equals 1.6ujq in this case, the phonon distribution func- 
tion remains well described by the antiadiabatic result 
which applies, strictly speaking, to T <^ ujq. Physically 
this means that the degree of undressing of the polaron 
is rather small for this particular value of T. A qualita- 
tive change in the profile of P(n) takes place upon going 
from R = 0.3 to R = 0.6, marking the progressive un- 
dressing of the polaron. Finally for large T, exemplified 
by R = 2 in the lower panel, the distribution function 
P(n) approaches the single-peak structure of Eq. ([57]) . 
which is plotted for comparison by the dotted line. As 
discussed in Sec. lIIlEl the phonon is too slow to track the 
individual tunneling events in this limit, hence it samples 
only the average displacement X(hd) = A/2. 

A useful perspective on the phononic state is pro- 
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FIG. 10: (Color online) The phonon distribution function 
P(n) = (\n}(n\), for U) /D = 0.0025, A = 6, R = Y/E p = 0.5, 
T — > 0, and different values of No- The detuning energy Ae d 
is held fixed at zero for all values of iVo, such that the elec- 
tronic level is half filled. NRG parameters: A = 2, N s = 4000, 
and N h = 1600. 
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FIG. 11: (Color online) Schematic description of a molecular 
bridge, where a single molecule is trapped between two leads. 

of n, it is clear that the low-energy state of the system 
can not be well represented using just one or two coher- 
ent states. This sets a stringent constraint on variational 
treatments of the crossover regime when A is large. 



vided by varying the reference charge No of the elec- 
tronic level. As discussed in Sec. [Til No can be elimi- 
nated from the Hamiltonian of Eq. (J]} by simply shift- 
ing the bosonic mode according to B = b — XN . If 
the level energy e d is simultaneously adjusted such that 
Ae d = e d — X 2 ujq(1 — 2Nq) is held fixed, then the varia- 
tion of No has no effect on the low-energy spectrum of 
the Hamiltonian. Nevertheless, the low-energy state of 
the system does vary with N , as can be seen from the 
case where the B boson occupies a pure coherent state. 
If B resides in the coherent state \z), then the b phonon 
occupies the shifted coherent state \z + XN ). In other 
terms, one can continuously vary the coherent state that 
b occupies by tuning No while holding Ae d fixed. This 
strategy can be used not only to expose a pure coherent 
state, but also to diagnose a superposition involving a 
small number of coherent states of the type often used in 
variational calculations. 

Figure [TU] displays the evolution of the phonon distri- 
bution function with increasing No, for A = 6 and fixed 
Ae d = 0. The ratio T/E p = 0.5 is tuned to the mid- 
dle of the crossover regime, such that the system is well 
removed from the extreme adiabatic and antiadiabatic 
limits. As can be seen, P(n) evolves in a rather com- 
plicated manner upon going from No = to No = 0.5. 
Initially, there are two rather broad and smooth humps 
for A^o = 0. Upon increasing No, the phonon distribution 
function narrows considerably, developing multiple sharp 
peaks by the time No = 0.5. Indeed, 3 clear maxima are 
visible at n — 0, 2, and 4 when No = 0.5, indicating 
that at least three independent coherent states signifi- 
cantly contribute to the low-energy state of the system. 
We have not succeeded in reproducing this multiple-peak 
structure for No = 0.5 (let alone the entire collection of 
distributions for the different values of No) by employ- 
ing a simple superposition of just a few coherent states. 
Given the sharpness of the peaks found at even values 



V. CONDUCTANCE OF A MOLECULAR 
BRIDGE 

Our discussion thus far has focused on thermodynamic 
quantities, which are difficult to measure for a single 
molecule. Of particular interest are transport properties 
of molecular devices of the type displayed schematically 
in Fig. II H where a single molecule is trapped between 
two leads. Continuing with spinless fermions, we model 
such a molecular bridge by the Hamiltonian 

n = E e kC { a c ka + E 4= E i dt ^ + H - c -} 

k,a k 

+ e d dU + oj tfb + Xoj (dU - N ) (6* + 6) , (49) 

where a = L,R labels the left/right lead. Here, ct cre- 
ates a conduction electron with momentum k in lead a, 
while t a is the tunneling matrix element between the elec- 
tronic level and the Wannier state closest to the molecule 
in lead a. For convenience we take the number of lattice 
sites N to be the same in both leads, although this con- 
dition can easily be relaxed. 

In equilibrium, the two-lead Hamiltonian of Eq. (|49)l 
is equivalent to the single-band model of Eq. (fT]). This 
stems from the fact that the localized level couples solely 
to the "bonding" combination 

while the "anti-bonding" combination 

is decoupled from the level. Inasmuch as impurity-related 
quantities are concerned, one can omit then the "anti- 
bonding" band altogether, to be left with the single-band 
Hamiltonian of Eq. (fT]) where t = y/t\+ t\. 
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FIG. 12: (Color online) The scattering phase shift S vs the 
detuning energy Ae d , for uio/D — 0.0025, A = 6, T — > 0, and 
three intermediate values of R — T/E p . Here V = Tl + Fr is 
the total tunneling rate from the level to the two leads. Solid 
lines depict the phase shift 6/ir, dashed lines show the dot 
occupancy rid, and the symbols show the number of displaced 
electrons, AiV. In accordance with the Friedel-Langreth sum 
rule, 8/n coincides with AiV, which deviates, however, from 
n d . NRG parameters: A = 2, N s = 4000, and N h = 1600. 



FIG. 13: (Color online) The zero-temperature conductance G 
vs the detuning energy Ae d , plotted for the same model and 
NRG parameters as in Fig. 1121 Here Go is a geometric factor 
specified in Eq. (|53[) . encoding the asymmetry of the tunnel 
junction. The characteristic shoulder in rid translates to a 
similar shoulder in G/Go, providing a distinct experimental 
fingerprint of phonon-assisted tunneling when the electron- 
phonon coupling is strong. 



The restriction to the "bonding" band is no longer 
complete when a finite bias is applied between the leads, 
since the nonequilibrium boundary condition applies to 
the lead electrons rather than their "bonding" and "anti- 
bonding" combinations. Thus, one can no longer settle 
with the single-band Hamiltonian of Eq. (JTJ for a biased 
junction. Fortunately, this complication can be circum- 
vented in linear response, when the conductance can be 
expressed in terms of equilibrium response functions of 
the system. In particular, the zero-temperature conduc- 
tance, G, is determined by the scattering phase shift of 
the "bonding" electrons, as follows from the Fermi-liquid 
ground state of the system. The zero-temperature con- 
ductance thus takes the standard form 



G = G sin 2 (S), 



where 



Gn — 



4r L r R 



h (Ti 



(52) 



(53) 



with r a = Trpot^ is a geometric factor encoding the de- 
gree of asymmetry in the coupling to the two leads and 
S is the scattering phase shift. Generalizing the Friedel- 
Langreth sum rule^ to the present setting, the phase 
shift S is given by 



6 = TrAiV, 



(54) 



where AiV is the number of displaced electrons. The lat- 
ter quantity comprises two contributions: (i) occupancy 
of the level n^, and (ii) the change in occupancy of the 
"bonding" band inflicted by the coupling to the level. 



In the wide-band limit only the former contribution is 
left, resulting in 8 = nrid- However, as we show be- 
low, AiV deviates from rid when the conduction-electron 
bandwidth is finite. 

Figure [T2l shows the scattering phase shift 5/ir versus 
Ae d , calculated directly from the NRG spectra using the 
standard prescription^ of Eq. (jBlip . For comparison, 
we also plot the level occupancy rid and the number of 
displaced electrons AiV, which is readily computed in the 
NRG since the total electronic occupancy of the system 
is a good quantum number that is kept track of in the 
course of the iterative procedure. The electron-phonon 
coupling is set equal to A = 6, while r = Tl + r# is 
tuned to three intermediate values of R = F/E p where 
rid shows a pronounced shoulder. 

While S/tt and AN are indistinguishable to within nu- 
merical accuracy, there is a small but consistent deviation 
from rid, particularly along the characteristic shoulder in 
rid- The difference between AA^ and rid appears to some- 
what increase upon going from R = 0.2 to 0.4, and dis- 
appears for Ae d —> (not shown). This latter behavior 
is to be expected since S/tt, AiV, and rid are all pinned 
to 1/2 by particle-hole symmetry when Ae d = 0. 

The corresponding zero-temperature conductance is 
shown in turn in Fig. [13] as a function of the detuning 
energy Ae d . Note that the conductance is an even func- 
tion of Ae d , as follows from the particle-hole transfor- 
mation of Eq. ©. Hence only the regime Ae d > is 
displayed. The characteristic shoulder in rid translates 
to a similar shoulder in G/Gq, which is further pushed 
down in magnitude due to the quadratic dependence on 
sin(<5). Experimentally, Ae d can be tuned using a suitable 
gate voltage, giving rise to characteristic shoulders in the 
conductance as a function of gate voltage. These shoul- 



14 



ders, which are quite unique in the context of tunneling 
through confined nanostructures, provide a distinct ex- 
perimental fingerprint of phonon-assisted tunneling when 
the electron-phonon coupling is strong. 

VI. SUMMARY 

Focusing on strong electron-phonon interactions, 1 <C 
A, we presented a comprehensive study of the crossover 
from the antiadiabatic to the adiabatic regime of phonon- 
assisted tunneling in the framework of a minimal model 
for molecular devices: a resonant level coupled by dis- 
placement to a single localized vibrational mode. Our 
main findings are as follows. 

1. In contrast to common lore, the crossover from the 
polaronic physics of the antiadiabatic limit to the 
perturbative physics of the adiabatic regime is gov- 
erned primarily by the polaronic shift E p — \ 2 u>o 
rather than the phonon frequency ujq. In particu- 
lar, the perturbative adiabatic limit is approached 
only as the bare hopping rate T exceeds the po- 
laronic shift, leaving an extended window of cou- 
plings where T well exceeds the phonon frequency 
and yet the physics is basically that of the antiadi- 
abatic regime. 

2. Throughout the traditional and the extended antia- 
diabatic regimes, the effective low-energy Hamilto- 
nian at energies below luq is the purely fermionic 
IRLM, which depends at resonance on two param- 
eters only: t c g and U c g. The effective tunneling 
amplitude i c ff obeys the empirical scaling form of 
Eq. (14*5)) . at least up to values of T several times 
larger than ujq, while C/ e ff is well approximated by 
Eq. (|2"6")) for much of the extended antiadiabatic 
regime. 

3. Although r c ff varies by many orders of magnitude 
as a function of A and R = T/E v , it is rather well 
described for all parameter regimes by the empirical 
formula of Eq. (|4T|) . which depends on a single scal- 
ing function F(R). Our proposal for T(R) is given 
by the solid line in Fig. [5j although this choice can 
possibly by further optimized. 

4. Charging properties are governed by two distinct 
mechanisms at the extended antiadiabatic and into 
the crossover regime. At small detuning, Ae d ^ 
r c g, the level occupancy rid follows a strongly 
renormalized noninteracting form with T — > T e f{. 
By contrast, a simple perturbative expansion in t 
applies for r c ff <C Ae d , giving rise to a characteris- 
tic shoulder in rid and in the low-temperature con- 
ductance of a molecular junction as a function of 
gate voltage. 

Our scenario for the crossover from the antiadiabatic to 
the adiabatic regime is summarized in Fig. [TJ where we 



have merged the traditional and the extended antiadia- 
batic regimes on the basis the two are qualitatively the 
same. 

This study was devoted to thermodynamic properties, 
for which a rather complete picture was provided. How- 
ever, the question of real-time dynamics remains largely 
open. Particularly, does r p = 1/E p show up as a new 
time scale in the dynamics? As could be anticipated, 
r p plays no role deep in the adiabatic regime, when 
E p <C T. Indeed, as recently shown for different scenarios 
of quench and driven dynamics, 35 only three time scales 
are involved in this limit: the dwell time r^wcii = 1/r, 
the period of oscillations r osc = 27r/oj c ff with the soft- 
ened frequency of Eq. ([2~Tj) , and the phonon damping time 
Tdamp = 7r r 2 /(A 2 a;o), extracted from the imaginary part 
of Eq. ([20"|) . Other perturbative2S and numerical^— 
studies of real-time dynamics have focused primarily on 
the build up of the current in a biased two-lead setting, 
revealing rich behavior. Most notably, a significant de- 
pendence on the initial state was reported in Ref. up 
to times far exceeding r^weii, suggesting the existence of 
another, much longer time scale whose origin is not quite 
clear. We note in passing that the parameter set referred 
to as adiabatic in the latter study corresponds in our 
notation to A = 8 and R = 1/8, which actually falls in 
the extended antiadiabatic regime. This may suggest the 
possible relevance at long times of a strongly renormal- 
ized tunneling rate akin to T ff . Clearly the understand- 
ing of real-time dynamics in molecular devices is a timely 
and challenging task that deserves further investigation. 
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Appendix A: Schrieffer- Wolff-type transformation 

for T < cjo 

In this Appendix, we describe the Schrieffer- Wolff-type 
transformation that maps the Hamiltonian of Eq. ([22]) 
onto the IRLM Hamiltonian of Eq. with the cou- 

pling constants t c g and U c g specified in Eqs. (1231) and 
(|26l) . respectively. As emphasized in the main text, the 
mapping applies to T < wq, when a single-step elimina- 
tion of all excitation energies exceeding ujq is sufficient. 
The mapping is further restricted to temperatures be- 
low ujq, when the (transformed) phonon is frozen in its 
unperturbed ground-state configuration b^b = 0. 

Following Schrieffer and Wolff,— we seek a canonical 
transformation 

H s = e s H'e- s (Al) 
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with a suitable anti-Hermitian operator S such that the 
low-energy subspacc is decoupled to order t 2 from all ex- 
cited phononic states and all single-particle conduction- 
electron excitations with energy |e&| > ujq (whether par- 
ticle or hole). To this end, we introduce two complemen- 
tary projection operators, P and Q = 1 — P, where P 
projects onto the low-energy subspace with Wb = and 
no single-particle conduction-electron excitations whose 
energy exceeds u>o- In other terms, all conduction- 
electron modes with energy e k < — uj (ujq < e k ) are 
strictly occupied (unoccupied) within the subspace de- 
fined by P. The Hamiltonian T~L' is then divided into an 
"unperturbed" part H' Q and a "perturbation" H[, where 



resulting in 



H' = PH'P + QH'Q 



and 



n[ = PU'Q + QH'P. 
Explicitly, T-L' takes the form 

restric 

+ ^E{ eA(fct ^ tc fc + H - c > 



with 



U o = J2 ekC l c k 



(A2) 



(A3) 



P 



(A4) 



(A5) 



while H[ is given by H[ =H[P + H.c. with 



Cfe 



n>0 fc 



clcc 



hole 



|n)(0| 
|0><0] 



and 



(±xy 



(A6) 



(A7) 



Here the symbols Y^k^ an d I]J°' C come to indicate that 
the summations over k and q are restricted to momenta 
such that 6^ > loq and t q < — loq, respectively. 
Using the formal expansion 

H s =n / Q +H' 1 + [S,H' } + [S,H' 1 }+^[S,[S,n' }}+... (A8) 

and anticipating that S is proportional to t at leading 
order (as will shortly be seen), one can group the dif- 
ferent terms in Eq. (IA8|) according to powers in t. The 
requirement that no coupling is left to order t 2 between 
the excited and low-energy subspaces is satisfied by de- 
manding that 



(A9) 



n s = n' + -[s,n' 1 ] + o(t 3 ) 



Equation (|A9|) has the formal solution 



1 



s = -—n' 1 + o(t 3 ), 



(A10) 



(All) 



where C is the Liouville operator defined by C' O = 
[O,%o]. Since is proportional to t, it is clear that S 
has a leading linear dependence on t, as we have assumed. 
The anti-Hermitian operator S does contain, however, 
additional higher order terms in t, which stem from the 
fact that Hq (and thus C' ) includes components linear 
in t [see Eq. (|A4[) ]. Denoting the linear-order compo- 
nent of S by the latter is computed by substitut- 
ing H' Q — > Ho, corresponding to setting C' — > Co with 
CoO = [d,Ho] in Eq. ([ATT]) . Using H[ of Eq. jA6|), this 
yields the explicit expression 



= SP- H.c. 



with 



S = 



V n>0 k 



c{d- 



u- 



Ck - 



t 



N 



€k + nujQ 

clcc -. hole -j 



nuo - £fc 



(A12) 



|n)<0| 



|0)(0|. (A13) 



To obtain Hs up to second order in t, it suffices to 
replace S in Eq. dATO)) with of Eqs. (IAT21) and (TA131) . 
Carrying out the commutator in Eq. (IA10[) . projecting 
the result onto the 6 T 6 = subspace, and restricting the 
band to \e k \ < uj , one arrives to order t 2 at the effective 
low-energy Hamiltonian of Eq. (pM)) with the coupling 
constants specified in Eqs. (|25|) and (|26l) . 

Two comments should be made about the derivation 
of the effective Hamiltonian of Eq. (|24l) . First, we have 
explicitly assumed a particle-hole symmetric band. Sec- 
ond, we have neglected Ck in the denominators on the 
first line of Eq. (| A13|) . on the premise that \ek\ is small 
as compared to the new ultraviolet cutoff energy ujq ■ 



Appendix B: Extracting the couplings of the 
effective IRLM Hamiltonian 



As discussed in the main text, the Hamiltonian of 
Eq. ([T]) with 1 <C A and T E p can be described at 
energies below uio by an effective IRLM. For Ae d = 0, 
the case considered hereafter, particle-hole symmetry re- 
stricts the number model parameters in the IRLM to 
two: the tunneling amplitude i e ff and the local Coulomb 
repulsion U c g. In this Appendix, we describe in detail 
how these model parameters can be extracted using the 



16 



Energy level 


Quantum numbers 


Hopping matrix 


(A = 2) 


(Qd,Qc) 


element to 


u 


(-1,0) 




(+1,0) 


— 




(-1,-2) 


— 


n 4Qi fi 


(+1,-2) 
C-1,+2) 






(+l,+2) 




0.9832 


(-1,0) 

(+1,0) 





TABLE I: Finite-size NRG spectrum at the noninteracting 
free-impurity fixed point t = U — 0, for A = 2 and odd iter- 
ation number N. The left column specifies the dimensionless 
NRG energies. The corresponding quantum numbers, Qd and 
Qc are listed in the central column, while the right-hand-side 
column indicates degenerate eigenstates that are connect to 
the designated state with a nonzero hopping matrix element 
once a finite t is switched on. 



NRG. The analysis relies on certain characteristics, ex- 
posed below, of the finite-size spectrum of the IRLM near 
the free-impurity fixed point, t e & — U c g = 0. 

In our analysis we shall use the IRLM in its following 
representation 



k V k k 

+ |_( d t d _ 1/2) £ :c t Cfc , : 



H.c.} 



(Bl) 



k,k' 



where is the number of k points that are summed 
over. This form of the Hamiltonian differs from that of 
Eq. (|24[) in the normalization factors multiplying t and 
U. In Eq. (|24l) . the summations over the momenta are 
restricted to N' < N values of k, hence the conversion 
between (t e s,U e g) and (t,U) reads as 



Energy level 


Quantum numbers 


Hopping matrix 


(A = 2) 


(Qd, Qc) 


element to 


n 

yj 


(-1,-1) 
(+1,-1) 


— 


(-1,+1) 


(+1,-1) 




(+1,+1) 






(-1,-3) 


— E 




(+1,-3) 


(-i,-i) 




(-1,-1) 


(+1,-3) 


0.9723 


(+1,-1) 
(-1,+1) 


— 




(+1,+1) 


(-1.+3) 




(-l,+3) 


(+1,+1) 




(+1.+3) 


— 


1.9446 


(-1,-1) 
(+1,-1) 
(-1,+1) 


(-1,+1) 
(+1,-1) 




(+1,+1) 





TABLE II: Same as Table [J for even iteration number N. 



1. Finite-size spectrum of the IRLM near the 
free-impurity fixed point 

a. Fixed-point spectrum for t = U = 

Consider first the noninteracting free-impurity fixed 
point of the IRLM with £ = {7 = 0. Since the impurity 
level is decoupled from the band, the fixed-point spec- 
trum is simply that of a free symmetric band, with an 
extra degeneracy of two due to the impurity level which 
can be either empty or full. Thus, the ground state is 
doubly degenerate for odd NRG iterations N, and four- 
fold degenerate for even iterations. The corresponding 
eigenstates are conveniently labeled by a pair of numbers 



Q d = 2tfd - 1 



(B4) 



In' tt tt n' 



(B2) 



For the linear dispersion considered here, the ratio N'/N 
equals D c g/D, where D e g ~ ujq is the effective band- 
width in the Hamiltonian of Eq. (|24l) and D is the full 
bandwidth pertaining to the original electron-phonon 
Hamiltonian of Eq. ([TJ. Thus, the conversion between 
the corresponding coupling constants becomes 



t = t, 



cir 



D 



err 



D ' 



U = U, 



D. 



err 



eff- 



D ' 



(B3) 



which is naturally accounted for, as we shall see, in the 
NRG. 



and 



Q c = ^[2/t/„-l] : 



(B5) 



n=0 



which serve as good quantum numbers for zero tunnel- 
ing (only their sum Qd + Q c is conserved for nonzero t) . 
The two degenerate ground states for odd N correspond 
to {Qd,Qc) = (±1,0), while the four degenerate ground 
states for even N are labeled by (Qd,Qc) — (±1,±1)- 
The extra two-fold degeneracy for even N stems from 
the presence of a conduction-electron mode that lies ex- 
actly at the Fermi level. The fixed-point spectra for odd 
and even iterations are listed in Tables Q] and |H] up to the 
second excitation energy. 
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FIG. 14: (Color online) Finite-size spectrum of the IRLM near 
the free-impurity fixed point, for U/D = 0.01, t/D = 10 -9 , 
A = 2, and even iterations TV. A nonzero U lifts the four- 
fold and eight-fold degeneracies of the ground state and the 
first excitation energy, respectively, splitting them into dis- 
tinct doublets (middle two panels). An infinitesimal t further 
splits some of the doublets. In particular, the ground-state 
doublet is split according to r y(t/D)A N< - 1+a ^ 4 ' (right panel, 
green dots), where 7 equals 0.835 and a is given by Eq. (|B7|) . 
The elementary particle, hole, and particle-hole excitations 
in the Qd = 1 sector are marked by rj+, r/_, and 2r]o in the 
middle panels (see text for details). 



b. Fixed-point spectrum for t — and U 7^ 

A small but finite U lifts certain degeneracies of the 
U = spectrum while maintaining the same pair of quan- 
tum numbers (Qd,Qc)- In particular, the ground-state 
quartet for even TV is split into two doublets, each com- 
posed of two particle-hole-symmetric states. The order of 
doublets discloses the sign of U. When U > 0, the states 
labeled by (+1, -1) and (-1, +1) [(+1, +1) and (-1, -1)] 
form the ground-state [excited] doublet, while the order 
is reversed for U < 0. The magnitude of U can be de- 
duced in turn from a standard phase-shift analysis^ (see 
below), which requires identification of the elementary 
particle, hole, and particle-hole excitations. For U > 
and even TV, these are given by the excitation energies rj + , 
?7_, and 2r] of the Qd — 1 sector, depicted in the middle 
panels of Fig.[l4j In the Qd = — 1 sector the elementary 
particle and hole excitations are interchanged, reflecting 
a reversal in sign of the scattering potential experienced 
by the band electrons when the level is empty. 



c. Spectrum in vicinity of the free-impurity fixed point 

In contrast to U, which is a marginal perturbation, 
the tunneling term t is a relevant one, driving the sys- 
tem away from the free-impurity fixed point to a new 
strong-coupling fixed point with 7r/2 phase shift of the 
scattered electrons. The new fixed-point spectrum dif- 
fers substantially from the t — one. However, it takes 
these differences a while to develop in the course of the 
NRG iterations. Starting with an infinitesimal tunneling 



amplitude, the finite-size spectrum displays only minis- 
cule deviations from the t — one over many iterations. 
This behavior persists as long as the renormalized tun- 
neling amplitude remains small. We focus hereafter on 
U > (the case relevant to phonon-assisted tunneling) 
and on this portion of the NRG level flow. 

The main effect of t at these iterations is to lift some 
of the remaining degeneracies of the t — fixed-point 
spectrum. For even TV, certain levels are split linearly in 
t, most notably the ground-state doublet. The splittings 
for odd TV are of higher order in t, reflecting the absence 
of a direct tunneling matrix element between degenerate 
eigenstates of the t = spectrum. (A detailed analysis of 
these matrix elements is presented in the right-hand-side 
columns of Tables HI and HT1) As we next describe, splitting 
of the ground-state doublet for even TV is proportional to 
t{Dfi) = t(D]y)/D]y, where t(D^) is the renormalized 
tunneling amplitude at energy Djq. 

According to a perturbative renormalization-group 
(RG) analysis of the IRLM^ the dimensionless tunnel- 
ing amplitude t — t/D obeys the RG equation 



dt 
~dl 



1 



(1 + a)t 



with 



4 (2 
a = —5xj - -<5tj 



(B6) 



(B7) 



Here, 5\j — arctan(7rp C//2) is the phase shift associ- 
ated with U in the absence of tunneling, po = l/(2-D) 
is the conduction-electron density of states, and I equals 
hi(D/D') with D' the running bandwidth. Note that 
Eq. (|B6|) is perturbative in t but includes all orders in 
U. It is supplemented in principle by a second RG equa- 
tion describing the renormalization of yet the latter 
equation can be neglected when t is small. Solution of 
Eq. (|B6I) yields the renormalized tunneling amplitude at 
energy D': 



t(D')=t(D) 



(l+ Q )/2 



(B8) 



In the context of the NRG, D/D' is replaced at iteration 
TV with Dq/D n = A N / 2 , resulting in 



t(D N ) = A"^")/ 4 ( 



(B9) 



We have found empirically that the ground-state split- 
ting for even TV is proportional to t(D]\[), provided the 
latter coupling is not too large. Here, by ground-state 
splitting we refer to the difference in energy between 
the ground state and the first excited state of the Q — 
Qd + Qc = sector [the latter state is no longer the low- 
est global excitationiS if t(D N ) > U/D N ]. The accuracy 
of our statement is demonstrated in the right-most panel 
of Fig. [HI where the ground-state splitting (red line) is 
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compared for U/D = 0.01 to jt(D N ) with 7 = 0.835 
(green dots). Excellent agreement is obtained. A similar- 
degree of accuracy extends to all values of U > 0, pro- 
vided t(Dff) < 0.1. For larger values of t(Djv), the linear 
relation between t(D^) and the ground-state splitting 
gradually breaks down due to the departure from weak 
coupling. As for 7, its value depends on both U and A. 
For A = 2, the discretization parameter used throughout 
this work, 7 grows monotonically from 0.835 to 0.94 in 
going from U = to U -> 00. For < U/D < 0.3, 
the regime of relevance to our discussion, 7 changes by 
no more than 4%, allowing one to use the single figure 
7 = 0.835 in order to extract t e g. 



2. Extracting the couplings i e g and [7 e g 

Our analysis thus far has provided us with a thorough 
understanding of the finite-size spectrum of the IRLM 
near the free-impurity fixed point. Next we specify how 
one can exploit this knowledge to extract the coupling 
constants t e g and U c g that enter the effective IRLM of 
Eq. 424]). 

As anticipated in Sec. MI CI the finite-size spectrum of 
the Hamiltonian of Eq. (fTJ) is found to be well described 
for r c ff <C Dn *C loq by the weak-coupling spectrum 
of the IRLM with U > 0. The sign of U is exposed 
from the nondegenerate ground state for even iterations 
N, which belongs to the Q — sector4i To extract the 
model parameters of the effective IRLM Hamiltonian we 
have implemented the following procedure. First, a par- 
ticular iteration number N = 22 was chosen such that 
Dn = 0.146a;o. The dimensionless tunneling amplitude 
t(Dff) — AE/j was next extracted from the energy 
splitting AE between the ground state and the first ex- 
cited state of the Q — sector. The value of 7 was 
set equal to 7 = 0.835, in accordance with the discus- 
sion above. Lastly, the dimensionless Coulomb repulsion 
U(Dn) = U{Dm)/Dm was extracted from the elemen- 
tary particle, hole, and particle- hole excitations 77+, rj-, 
and 2?7o according to the standard prescription^ 



with 



6 = 



V+ 



7T 

2' 



(Bll) 



The procedure outlined above provided us with esti- 
mates for t and U at the energy scale D e ff = Dn=22- 
The accuracy of the couplings so obtained can be sum- 
marized as follows. If less than 0.1, the dimensionless 
tunneling amplitude t is accurate to within about 4%, 
provided U is simultaneously smaller than 0.3. This cri- 
terion was met for all points displayed in Fig. [5] Accuracy 
of the dimensionless Coulomb repulsion U is controlled 
in turn by the ratio t/U. Indeed, Eqs. (IBlOj) and (IB11|) 
are precise for t = 0, acquire a small correction propor- 
tional to t/U when t is small, and break down as soon as 
t approaches U. All points displayed in Fig. [7] fall in the 
range t/U < 2 x 10~ 3 , corresponding to an accuracy of 
order 1% for U. 

Finally, the conversion from the dimensionless cou- 
plings t and U to t e g and U e g proceeds as follows. Using 
the notation of Eq. (|B3|I , one has that 



/ii=/ ^^ =jD W^ 



(B12) 



and 



C/eff 



D 



D, 



-U = DU. 



c rr 



(B13) 



Since D e g /D is replaced at the Ath iteration of the NRG 
with Dn/D — A _7V / 2 , we arrive at 



D 



D 



(B14) 



U 



■ tan((5), 



(B10) 
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